Joint Galaxy-Lensing Observables and the Dark Energy 
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Deep multi-color galaxy surveys with photometric redshifts will provide a large number of two- 
point correlation observables: galaxy-galaxy angular correlations, galaxy-shear cross correlations, 
and shear-shear correlations between all redshifts. These observables can potentially enable a joint 
determination of the dark energy dependent evolution of the dark matter and distances as well as 
the relationship between galaxies and dark matter halos. With recent CMB determinations of the 
initial power spectrum, a measurement of the mass clustering at even a single redshift will constrain 
a well-specified combination of dark energy parameters in a flat universe; we provide convenient 
fitting formulae for such studies. The combination of galaxy-shear and galaxy-galaxy correlations 
can determine this amplitude at multiple redshifts. We illustrate this ability in a description of the 
galaxy clustering with 5 free functions of redshift which can be fitted from the data. The galaxy 
modeling is based on a mapping onto halos of the same abundance that models a flux-limited 
selection. In this context and under a flat geometry, a 4000 deg 2 galaxy-lensing survey can achieve a 
statistical precision of o(Q,ue) = 0.005 for the dark energy density, <t(wde) = 0.02 and u(w a ) = 0.17 
for its equation of state and evolution, evaluated at dark energy matter equality z w 0.4, as well as 
constraints on the 5 halo functions out to z = 1. More importantly, a joint analysis can make dark 
energy constraints robust against systematic errors in the shear-shear correlation and halo modeling. 



I. INTRODUCTION 

In the successful standard cosmological model where 
structure in the universe originates from Gaussian ran- 
dom density fluctuations in the initial conditions, all sta- 
tistical properties of cosmological structure observables 
depend on a single quantity: the linear power spectrum 
of mass fluctuations. The evolution of this spectrum de- 
pends on the properties of the dark energy in a precisely 
calculable way. The task of extracting dark energy in- 
formation from cosmological structures reduces to deter- 
mining the relationship between observables and the un- 
derlying linear mass power spectrum. 

Deep, multi-color photometric galaxy surveys measure 
two sets of cosmological observables: the angular dis- 
tribution of the galaxies and the weak lensing shear in- 
duced on their shapes. From these observables, three 
types of two-point correlations can be constructed: the 
angular correlations between the positions of the fore- 
ground galaxies, the shear-shear correlations between 
background galaxies, and the galaxy-shear cross corre- 
lations induced by the association of dark matter with 
foreground galaxies. 

Whereas these correlations have so far been analyzed 
separately and/or with data from different surveys, the 
joint analysis of these measurements will be feasible from 
forthcoming surveys. In this paper, we consider what can 
be learned from the combined two point correlations. We 
shall see that with the multitude of observables avail- 
able, prospects for the joint determination of the cosmol- 
ogy and the relationship between galaxies and mass are 
bright. 

Galaxy-shear cross correlations, also known as galaxy- 
galaxy lensing correlations, were first detected by Brain- 



erd et al. [[J, following earlier upper limits 0]. The 
Red-Sequence Cluster and VIRMOS-DESCART surveys 
(0, also 3) are examples of the current state of the art 
in deep surveys. Though shallower, the Sloan Digital Sky 
Survey (SDSS) complements these given the largerpop- 
ulation of foreground and background galaxies [a, S • 
These observations have been interpreted in terms of 
the dark matter distribution associated with the lensing 
galaxies and their environment under the halo model for 
galaxies USE!. 

Ongoing and future surveys such as the Canada France 
Hawaii Legacy survey, Pan-STARRS, LSST, and SNAP 
[ll| will produce multi-color catalogs of galaxies. Photo- 
metric redshifts can then be estimated for these galaxies, 
allowing for measurements with foreground galaxies ex- 
tending to z ~ 1 and background galaxies up to twice as 
far in multiple redshift bins. 

Photometric redshifts significantly augment the 
prospects for joint galaxy-lensing studies beyond the cur- 
rent state-of-the-art. The number of observable cross- 
correlation functions between the galaxies and the shear 
scales as the product of the number of lens and source 
redshift bins and can easily exceed the number of galaxy- 
galaxy clustering observables and independent shear- 
shear correlation observables. 

We consider prospects for constraining the evolution 
of the dark energy through such surveys. The combi- 
nation of galaxy-shear and galaxy-galaxy correlations is 
particularly fruitful in that they allow for a joint solution 
of the evolution of the matter and galaxy distributions. 
These determinations can be cross checked against those 
from the shear-shear correlation. The latter depend only 
on the mass power spectrum but are typically subject to 
more severe systematic uncertainties. 

We begin in fjn] with a description of the statistical 



2 



methods employed. These methods may be applied to 
any model of the two-point correlations. In the Appen- 
dices, we describe the parameterization of the linear mass 
spectrum in the standard cosmological model, including 
convenient fitting functions for the dark energy effects, 
and a generalization of the halo model for the associa- 
tion of galaxies with the mass. Our model for the lat- 
ter is based on recent developments in galaxy simula- 
tions which rely on matching the observed number den- 
sity of galaxies to the predicted number density of halos 
[l2|. Utilizing these parameterized statistics, we study 
prospects for constraining the properties of the dark en- 
ergy in a spatially flat universe in 3IIII We conclude in 

m 



II. STATISTICAL METHODS 

In this section, we describe the basic statistical ap- 
proach to the joint study of galaxy and lensing power 
spectra. These methods are independent of the specific 
parameterization of the power spectra described in the 
Appendix and employed in the following section. We 
begin with a brief review of the relationship between an- 
gular and three dimensional power spectra in S HI Al Wc 
relate these to the traditional galaxy-galaxy lensing ob- 
servables in 311 Bl In 311 Cl we describe how the sensitivity 
of power spectra to underlying parameters may be quan- 
tified. 



A. Power Spectra 

Under the assumptions of statistical isotropy and small 
angles, the two-point observables of a set of two dimen- 
sional scalar fields xi (n) , where n represents the direction 
on the sky, is given by their angular power spectra 

(x* (1)^(1')) - (2ir) 2 8(l-l')C?*> , (1) 

where 1 is the Fourier wavevector or multipole 

dH 



Xi(n) 



(2tt) 



(2) 



Here and throughout we use the same variable to repre- 
sent the field in angular and Fourier space. 

Suppose these angular fields are related to three di- 
mensional source fields Si(r; z) by a weighted projection 



Xi(n) = J dzWi(z)si(ri = hD A ;z), 



(3) 



where D A (z) is the angular diameter distance in comov- 
ing coordinates. We will use semicolons to denote argu- 
ments that will be suppressed where no confusion will 



The Limber approximation 0, then relates the 
two dimensional power spectra to the three dimensional 
power spectra as 



C 



dz^ r W i {z)W j {z)P'*'i{k = l/D A ;z) , (4) 

A 



where H(z) = a 1 da/dt is the Hubble parameter and 

( S *(k) Sj (k')) = (2nfS(k-k')P^(k), (5) 

defines the three dimensional source power spectrum. 

For our purposes, the two dimensional fields will be 
the "lens" galaxy number density fluctuations and the 
electric or e component of the weak lensing shear field 
measured with "source" galaxies. The lens and source 
galaxies of a given survey may be divided into bins ac- 
cording to redshift, luminosity, color or other criteria. 
We will here concentrate on redshift binning. However 
for notational convenience we will suppress the subscripts 
i denoting the bins for the rest of this section. 

For the galaxy fluctuations, the source field is the three 
dimensional number density ny(r; z) or rather its fluctu- 
ations 



s(r;z) = S g 



5n v 

fly 



(6) 



and the weight for the angular fluctuation field g(h) is 
the normalized redshift distribution function 



W g (z) 



D\ n v 



where the normalization factor 

( j D A- 



(7) 



(8) 



is the angular number density in sr . Note that the 
weights are normalized so that J W g [z)dz — 1. 

For the weak lensing shear, the underlying scalar field 
is the electric component of the shear field x(fi) = e(n) 
which manifests itself as a shearing of background galaxy 
images according to the complex shear 



7i (n) ± «7 2 (n) 



d 2 l 



(9) 



where 4>i is the azimuthal angle of the Fourier vector with 
respect to the ei axis. The e field itself is a projection of 
the mass density fluctuation 



s(r; z) = S m = , 



(10) 



and hence is equal to the convergence K(n). 

The weights are given by the efficiency for lensing a 
population of source galaxies 

W e (z) = -fi m j J dz jy—W g {z ) , (11) 
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where the angular diameter distance to the lens is Dol — 
Da(D), to the source Dqs — Da(D') and between the 
lens and the source D^s = Da(D' — D). Here D(z) is the 
comoving coordinate distance and note that D = Da in a 
flat universe. The distribution of source galaxies W g {z') 
need not be the same as for the (lens) galaxies above. 
Furthermore W g , the normalized redshift distribution, is 
the direct observable so that the efficiency W e for a known 
W g may be used to probe cosmology. 

The complete two point statistics of the shear and 
galaxy correlations are thus specified by a choice of cos- 
mology and a description of the underlying three dimen- 



sional power spectra P' 



■S m 8„ 



P 



and P 



s a s„ 



as a func- 



tion of wavenumber k and redshift z. The latter two 
will depend not only on cosmology but also on galaxy 
properties. 



B. Cross Correlation Functions 

The galaxy-shear cross power spectra are related to 
the more familiar cross correlation functions through the 
Fourier transform relations Eqn. , © 



(7i(A)<5 9 (n')> 



(72(A)<5 9 (n')> 



d 2 l 

Idl 
2^ 
Idl 



Cf cos(20 z )e ll - (A - fl ' ) 



Cf e cos(2</>) J 2 (W), 



2tt 



Cf Bin(20) J 2 {W). (12) 



where n — n' = (9, cf>) is the angular separation vector 
with magnitude 9 and azimuthal angle <f> with respect to 
the e x axis. We have used the identity 

oo 

e ii-(n-n') = Jl (iQ) + 2 i m J m {W)cos[m{(f>i -</>)] . (13) 

m— 1 

Note that the correlation functions depend on both the 
magnitude of the separation vector 9 and the azimuthal 
angle <f>. Despite this complication, Eqn. (|12|l can be 
straightforwardly used to generalize a maximum likeli- 
hood estimator for the shear-shear angular power spec- 
trum (e.g. [H E3). 

Although the galaxy-shear cross power spectrum is 
thus a direct observable, observations to date have fo- 
cused on the tangential shear component around galaxies 
due mainly to systematic effects in the shear measure- 
ment. It is therefore useful to relate the two approaches. 

We can express the tangential shear about a galaxy at 
the origin as 



7 T (n) = -71 (n) cos(2</>) - 72(A) sin(20) 



(14) 



The angular correlation function then becomes a function 
of 9 alone and is given by 



(7r(0)>haio = (M0)S g (O)) 



(15) 



Under the ergodic assumption, this quantity can be rein- 
terpreted as the average tangential shear around lens 
galaxies in the sample volume. 

For a narrow redshift distribution of source and lens 
galaxies, the tangential shear directly probes the galaxy- 
mass power spectrum at the lens redshift. Substitution 
of the Limber equation (@J in Eqn. Ijl5|l gives 



(7t(t^-)) 



D, 



halo 



OL 



dk 

T 



Ho 
k 



k 3 P 



■5 a 6„ 



2tt 2 



(16) 



-MkR), 



where Po/Hq = 832Sl m ft,Af Q pc~ 2 and R is the distance 
transverse to the line of sight. The critical surface density 
is given by 

! _ ^G D OL D LS 3g 2 o DqlDls , 17 , 

2-i cr ~ — 9 r; — ~n — iim ■ l-'-'J 

c o.lP>os 2 p a L L>os 

For a single lens galaxy with an azimuthally symmetric 
density profile p(R, D), these quantities are related to the 
projected or surface mass density 



E(fl) = J dDp(R, D) , 



(18) 



through AE(iZ) = S(i?) — S(_R), where the average is 
over transverse distances interior to R. Note that all 
distances are in comoving coordinates. 

The tangential shear technique throws away informa- 
tion by combining the two components of the shear before 
averaging. It contains the majority of the signal when 
considering a spherically symmetric density distribution 
about a galaxy. At large radii, the usual approach is to 
attempt a reconstruction of the scalar e(n) (= k(A), the 
convergence in the weak lensing limit) out of both com- 
ponents (e.g. This reconstructed field then acts as 
a template density map for the correlation 



(e(0)6 g (0)) 



In 



(2*) aG ' e 



(19) 



Alternately, one can go the other way and define a tem- 
plate shear [l^ 



7f(n)±i7f(n) = 



d 2 l 

(27T) 



and construct the correlation 



2^(l)e ±2i *e fl - fi , (20) 



Idl 



(7?(%i(o)+7!(% 2 (o)) = J —crMie) 

= (e(6)6 g (0)). (21) 

Formally the two techniques construct the same correla- 
tion function and preference for one versus the other is a 
matter of considering errors in the reconstruction. 
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C. Fisher Matrix 

Given angular power spectra that are defined by a set 
of cosmological and galaxy parameters p a , forecasts on 
how well such parameters can be extracted from the data 
is an exercise in error propagation. 

In general, the observed two points statistics of the 
angular fields C; will receive a contribution from noise 
sources which we will assume to be statistically isotropic 



(22) 



We will further assume that the noise contributions for 
the galaxy and shear fields arise from uncorrelated shape 
and shot noise 



N 



7 2 

I rms 



Nf a9b = S ab ^- 

nAa 



N u 9a 



0. 



(23) 



where 7 rms is the rms shear in each component arising 
from the intrinsic ellipticity of the galaxies and measure- 
ment noise. 

The shear fields are expected to be nearly Gaussian 
with respect to power spectrum errors for I < 10 3 due 
to linearity and projection |19| as borne out in simu- 
lations |20|. For the galaxy fields, the transition scale 
is at somewhat lower I depending on the width of the 
projection (e.g. 0, H^)- However for both, the noise 
contributions at I 3> 10 3 dominate the errors and mask 
the non-Gaussianity of the underlying fields. Under the 
assumption of Gaussian statistics for the fields, the infor- 
mation contained in the power spectra can be quantified 
by the Fisher matrix 



(24) 



where the sum is over bands of width Al in the power 
spectrum and / s k y is the amount of sky covered by the 
survey. The rough / s k y scaling is valid for contiguous 
regions with comparable extent in each of the two angular 
directions. 

Here we have suppressed the indices in a matrix 
notation and 



ID,, 



D 



dc; 



dp 



(25) 



The inverse Fisher matrix approximates the covariancc 
matrix of the parameters C p ss (F _1 ). 

One can also break this information down into subsets 



F; 



sub 



a/3 



f sky J2(2l + l)Al 



j-^XiXj r^subl-l n . 

U la IM \ U , 



{ij)(mn) 



(26) 



where sum over ij and run pairs can run over a restricted 
subset of the observables. The covariance matrix of the 
subsetted power spectra is given by 



^£isub | ij,mn ^jXiX m Q x . 



c 



(27) 



In the limit that the sum is over all combinations, 
Eqn. H26f) returns Eqn. I|24|l . This subsetting also clar- 
ifies the role of the Gaussian assumption. Gaussianity 
implies a diagonal covariance matrix in I and reduces its 
form to the product of power spectra in Eqn. (|2*7Jl . Non- 
Gaussianity from non-linear structure formation corre- 
lates the power between high / bands but in a fairly sim- 
ple way: all bands share a common normalization whose 
variance is determined not by Gaussian statistics but by 
the sample variance near the non-linear scale. Under the 
halo model of Appendix B, this behavior arises because 
the shape reflects the shape of halo profiles whereas the 
amplitude reflects their abundance. This abundance fluc- 
tuates with the local mean density. 

Under the Limber approximation of Eqn. given 
Nl lens galaxy samples in disjoint redshift bins and N3 
source galaxies samples, there are Nl distinct galaxy 
spectra, N$(Ns + l)/2 shear spectra, and N^Ns galaxy- 
shear cross spectra. The potentially large number of cross 
spectra offers great opportunities for studies of galaxy 
evolution and cosmology. Consider then the Fisher ma- 
trix of cross spectra alone 



F ge — f 
r a0 ~ J'- 



sky 



^(2^ + l)AZ 



NlNs ara^i 

Mew) 



where 



QQagbQ^ 



<Cf> 



dp a 

(28) 
(29) 



Note that the total variance of the shear and galaxy fields 
contributes to the noise of the cross correlation as a type 
of sample variance. Furthermore in the low signal-to- 
noise regime, the covariance is dominated by the product 
of power spectra of the galaxy and shear fields not the 
sample variance of the signal. In this case, the Z-diagonal 
form of the Fisher matrix depends on the assumption 
of statistical isotropy and remains valid even when the 
shear and galaxy fields are strongly non-Gaussian. 

The framework described here is general and may be 
applied to any parameterized model for the underlying 
three dimensional power spectra P SmSm , p s 9 s m anc j pS g s fJ 
and any selection criteria for the galaxies. In Appendix A 
we describe the well-tested standard model for the linear 
mass spectrum and how it depends on the dark energy. 
In Appendix B we develop the halo model for the galaxy 
and cross spectra. Motivated by recent simulations which 
associate galaxies with substructure in dark matter ha- 
los, we utilize 5 free functions of redshift to describe the 
occupation of galaxies in dark matter halos. By discretiz- 
ing these functions into observed redshift bins we obtain 
a halo model parameterization with 5Nl parameters 
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Because even this multi-dimensional halo model may 
not be sufficiently realistic, we will attempt in the next 
section to separate cosmological information that does 
and does not depend on the details of the halo model. 
For instance in the large scale regime where galaxies clus- 
tering is nearly fully correlated with the mass, a mea- 
surement of the galaxy auto and cross power spectra is 
essentially a measurement of the mass power spectrum. 
Furthermore, the angular cross spectra as a function of 
source galaxy redshift scale with distance in a known way 
through the lcnsing efficiency for any choice of the under- 
lying three dimension power spectrum P S 9 S ™. 



III. DARK ENERGY PROSPECTS 

In this section, we study the prospects for dark energy 
constraints with galaxy-lensing power spectra. We begin 
in flIII Al by defining the fiducial survey and the galaxy 
selection. In Hill Bl we study the distance related or halo 
model independent information in the galaxy-shear cross 
correlation and in frill CI the joint constraints from all 
power spectra. 



A. Fiducial Survey 

For illustrative purposes, let us define a fiducial survey 
that is loosely based on a next generation ground based 
lensing survey. We take a source redshift distribution of 
the form 



W g - S cx z 2 e -(z/zw)2 



(30) 



with zw corresponding to a median redshift z me d;S = 1-5 
and an angular number density of t%a = 70 arcmin -2 . 
This corresponds roughly to a magnitude limit of / w 27. 
For the shape noise we take a shear rms per component 
of 7rms = 0.3 which reflects the intrinsic cllipticity and 
measurement errors of ground based observations (see 
e.g. Note that the noise variance scales as 7^ ms /nA 

in Eqn. We take a survey area of 4000 deg 2 . 

From the survey galaxies we choose a galaxy lens pop- 
ulation from the high luminosity tail. To balance sig- 
nal strength and halo model robustness per lens galaxy 
against lens abundance, we choose lens galaxies with an 
abundance corresponding to a mass threshold of M t h — 
1O 13 ' 5 /i _1 M0 in the fiducial model. Due to this tradeoff, 
the net signal-to-noise ratio is only weakly dependent on 
the threshold. Note that the mass threshold is not held 
fixed as halo and cosmological parameters are varied (see 
Appendix B) The fiducial mass threshold simply defines 
the redshift distribution and angular density of the lens 
galaxies. These are the quantities held fixed under the 
variations. 

For the redshift binning, wc typically choose between 
2-5 photometric redshift bins in the source distribu- 
tion. Further divisions do not substantially enhance con- 
straints on the dark energy given the broad efficiency 



~i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i 

sampling errors 

no sampling errors 




FIG. 1: Lensing efficiency ratio and 68% CL forecasts on a two 
parameter model for the dark energy (SIfje, u»de)- With two or 
more source galaxy populations, here Ng = 2, the ratios of galaxy- 
shear power spectra provides a measure of distances through the 
efficiency ratio. Here we take Nl = 10 lens galaxy populations 
and marginalize the amplitude of the galaxy-mass power spectrum 
in all bandpower and redshift bins. The inner contour shows the 
effect of neglecting the sampling errors and shows that sample and 
shot noise are comparable for a shear noise of 7 rm s = 0.3 and a 
source density of ua = 70 arcmin -2 . These parameters and a 4000 
deg 2 survey are assumed here and throughout. 



function [23]. These are taken to have equal extent in 
redshift out to 2z mec j ; 5 but with the last bin containing 
the remaining high redshift galaxies. For the lens galax- 
ies, we limit the populations to z < 1 since they will re- 
quire more accurate photometric redshifts. We typically 
take Nl = 10 lens galaxy bins reflecting a photometric 
redshift accuracy of Az — 0.1. The total lens popula- 
tion then has an angular number density of 0.026 gal 
arcmin -2 with a median redshift of z mo d;L = 0.7. 

Finally, we allow for uncertainties in the shape of the 
mass power spectrum and initial normalization by taking 
priors of a(hiS^) — a(n s ) — a{\nQ,} } h 2 ) — a{\a.Q, m h 2 ) — 
0.1 corresponding to a conservative interpretation of cur- 
rent constraints (see Appendix A and e.g. [12] )• With 
the fiducial lens survey, dark energy results depend only 
weakly on these prior assumptions. Note that we take 
no prior constraints on the dark energy parameters so 
that the projected constraints reflect only the potential 
of galaxy-lensing power spectra. 



B. Model Independent Constraints 

Comparing lensing observables for different source red- 
shifts has been proposed as a way of measuring distances 
with galaxy clusters [2^, l26l I27I |28| . Recently, wide field 
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weak lensing statistics have also been developed to ex- 
ploit this test [TBI |2^| . In the limit of infinitesimal lens 
galaxy redshift bins, i.e. in Eqn. (fTTJl 

W g , L -> 5(z - z L ) , (31) 

the ratio of the galaxy-shear cross correlation in multiple 
shear source bins depends only on the ratio of efficiencies. 
This in turn depends only on angular diameter distances 
and redshifts. In terms of the cross power spectra 

CP W e2 (z L )' 1 ' 

Hence the galaxy-shear cross correlation provides infor- 
mation on the dark energy which is immune to uncer- 
tainties in the modeling of the underlying galaxy-mass 
correlation. The ratio is also immune to sample variance 
in the signal, i.e. galaxy to galaxy variations in the under- 
lying mass correlation and hence the non-Gaussianity of 
the signal. In the Fisher matrix of Eqn. (|28() , if the noise 
terms Cf 9 and Cf e vanish, the covariance matrix becomes 
singular implying infinitesimal errors in distance param- 
eters. In fact the technique works for individual galaxies 
in principle. Unfortunately with realistic noise estimates, 
the tradeoff between model independence and sensitivity 
is severe. For example, a 10% change in the equation of 
state parameter u>de typically yields an ~ 0.1% change 
in the efficiency ratio and so measurement of the effect is 
only possible with large galaxy and shear surveys. Con- 
trast this with the several percent change in the absolute 
growth or shear amplitude implied by Eqn. (|A9I) and il- 
lustrated in Fig. IT2l 

Even for this statistic the halo model enters in three 
ways: by defining the strength of the signal, the sam- 
ple variance of the noise and the accuracy to which the 
redshifts of the galaxy lenses must be known. The sam- 
ple variance arises from contributions to the shear from 
structure along the line of sight not associated with the 
galaxy population and from the intrinsic clustering of 
galaxies. Furthermore, with finite-width redshift bins in 
the lens galaxy distributions, the efficiency ratios depend 
on the model for the evolution of the underlying power 
spectra across the bins ^BJ . Even with our 5 halo param- 
eter model for the evolution in each bin, the constraints 
are compromised. No constraints are possible with pho- 
tometric redshift bins of Az « 0.1 if we marginalize all 
5 parameters for a given bandpower in I. By combin- 
ing multiple bandpowers, one can recover dark energy 
information, but that amounts to utilizing information 
from the shape of the underlying correlation and again 
degrades the model independence of the effect. We will 
return to this type of constraint in the next section. 

To study the efficiency ratio test, we instead marginal- 
ize a single halo parameter per redshift and bandpower 
bin. This model is equivalent to marginalizing a con- 
stant amplitude or bias per bandpower in the underlying 
power spectrum. Hence it eliminates information from 
the shape but retains an assumption on the evolution of 



the galaxy-mass power spectrum. Results do not depend 
on the number of bands employed so long as the signal 
and noise dominated regimes are separated since the ad- 
ditional parameters are employed simply to remove shape 
information from measurements of different bands. We 
take the parameter to be the satellite normalization A s 
for convenience and have verified that the results do not 
depend on this choice. 

In Fig. ^ we show the constraints in the constant equa- 
tion of state, dark energy density (tode-Hde) plane for 
Nl = 10, TVs — 2 (see Appendix A for dark energy pa- 
rameter definitions). Note that when moving to a varying 
w, the errors on equation of state at the best constrained 
redshift u> p i vo t remain the same <r(w p i vo t) = 0"(iode), 
whereas those on SIde increase. We can also assess 
the effect of sample variance by artificially removing the 
terms in the covariance (|27() that are proportional to 
C"C ; ". This causes an improvement in the errors by 
~ \/2 indicating that the sample variance is compara- 
ble to shot variance. This near equality reflects the fact 
that these halos have a projected scale radius of order 1' 
(see Fig. Illll . On these scales, the assumed shape noise 
7rms I \/nAg ~ 0.02 is comparable to the cosmic shear. 
We include the information from all scales; in practice 
the signal to noise has converged well before our numer- 
ical maximum I = 30000. 

The discrepancy of our results with those of |29j is due 
to a factor of two error in the calculation of the signal 
(since corrected, astro-ph/0306046 v3), to a more real- 
istic model for halo profiles, and the inclusion of sam- 
pling errors. Comparison of the results also requires 
implementing their prior on IIde and their shape noise 
specifications. These results agree with an independent 
and concurrent study which also extended the model- 
independent techniques to shear-shear correlations [30l |. 

Raising the number of source (N$) and lens (N^) di- 
visions or choosing less rare galaxy tracers only slightly 
improves these constraints. Note that the latter entails 
employing less massive objects with higher number den- 
sity but weaker shear signal. Finally note that the tem- 
plate technique of |l5j | can be reexpressed as a measure- 
ment of the zero lag correlation of Eqn. I|21(l or a single 
bandpower 

7 - = / iB? cr ■ (33) 

Although the efficiency ratio constraints are relatively 
weak, they are fairly robust to both the halo model and 
non-Gaussianity of the correlation on small scales. More- 
over, they become more powerful when combined with 
complementary probes of shape and amplitude of the var- 
ious power spectra on large scales as we shall now see. 

C. Model Dependent Constraints 

Given the underlying halo model parameterization, 
cosmological and halo model parameters can be jointly 
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FIG. 2: Galaxy-shear power spectrum 68% CL constraints on a 
two parameter dark energy model (Orjg, hide) with power spec- 
trum information out to lp = 1000 (upper) and 3000 (lower). 
Galaxy-shear constraints (dashed lines) are complementary to the 
efficiency ratio test (this solid) and are assisted by the addition of 
galaxy-galaxy constraints (thick solid) which help determine the 
5N^ = 50 halo parameters. Joint constraint (shaded) is only 
weakly dependent on lp and hence non-Gaussian errors. Here the 
number of source redshift distributions N$ = 2. 



fit to the observable power spectra. Let us first con- 
sider constraints that are based on the galaxy-mass cross- 
correlation alone. In Fig. [21 (dashed lines) we show 
the cross power constraints for constant wde-^de and 
N L = 10, Ns — 2. We marginalize 5N L = 50 halo pa- 
rameters and vary the maximum I employed in Eqn. i|28|) 
[see Appendix B] . Let us focus on the opposite region to 
the previous section, that of large scales I < lp ~ 10 3 
where non-Gaussianity in the fields and inadequacies in 
the halo model are minimized. Under the halo model 
prescription, there is sufficient information in the power 
spectra to constrain a degenerate combination of wde 
and f2r>E- Note that the errors in the best constrained 
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FIG. 3: 68% CL constraints in a three parameter dark energy 
model («)pivoti w a, SIde) f° r shear-shear correlations only (light 
shaded ellipse) and all 2-point correlations (dark shaded ellipses) 
with Nl = 10 and Ng = 5 and lp = 3000. Also shown for com- 
parison are the efficiency ratio constraints (dashed lines) and joint 
galaxy-shear and galaxy-galaxy correlations (dotted lines). Even 
after marginalizing over hN^ = 50 halo model parameters, galaxy- 
shear with galaxy-galaxy power spectra have comparable constrain- 
ing power on the dark energy as shear-shear power spectra. Errors 
between iu p i V ot and w a are uncorrelated by definition. The pivot 
point z p j vot = 0.36 and is close to the epoch of dark energy domi- 
nation. 



direction are insensitive to the maximum / p and hence to 
uncertainties in the non-Gaussianity of power spectrum 
errors. They are also then less sensitive to inadequacies 
in the halo modeling that appear on small scales. 

The degeneracy line follows a line of constant linear 
power spectrum amplitude and distance at the typical 
redshift of the lenses. Breaking this degeneracy then 
depends both on the internal or external determination 
of parameters in the halo model and its overall validity. 
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since the J7de dimension is shown. We have chosen here 
Nl = 10 and N$ — 5. However shear-shear correlations 
are more susceptible to systematic errors since the cross 
power would null out systematics that are not common 
to both. They also require a modeling of the statistics 
and their covariance in the non-linear regime. Here we 
have taken lp — 3000. 



The combined galaxy-shear and galaxy-galaxy power 
spectra potentially have constraining power that is com- 
parable to the shear-shear spectra. The combination 
of all three reduces the errors by of order y/2 or more 
versus shear alone, leading to marginalized errors of 
o-(^de) = 0.005, ff(w pivot ) = 0.02 [a(w ) = 0.05], and 
&{wa) = 0.17. More importantly they provide important 
cross checks against shear systematics on the one hand 
and inadequate halo modeling on the other. 



FIG. 4: Example of combining external constraints. The 
Wpivot galaxy-lensing constraints can be transformed into 
other dark energy parameterization conventions (here Qde, 
wo, w a ) for comparison and joint studies. Dashed lines rep- 
resent the improvement to the 68% CL region due to the 
addition of projected CMB constraints for the Planck satel- 
lite which mainly constrain the angular diameter distance to 
recombination. 



Fortunately, the efficiency ratio information from small 
scales is complementary in direction. Furthermore, infor- 
mation from the other spectra can included. In partic- 
ular the galaxy-galaxy correlations are more sensitive to 
the assumptions of the halo model than the galaxy-mass 
correlation (see Fig. IT2")l . It may be used to cross check 
and calibrate the halo model parameters and potentially 
extend the parameter space as needed. 

In Fig. |21 we illustrate the utility of a joint analysis. 
Adding in galaxy-galaxy power spectra information out 
to the same lp substantially assists dark energy param- 
eter constraints by effectively acting as a prior and con- 
sistency check on the halo model parameters. The small 
interior ellipse illustrates the potential for simultaneous 
determinations of a constant wde and Ode by combining 
this information with that from the efficiency ratio. 

Just as galaxy-galaxy correlations provide a power- 
ful cross check on halo model parameter determinations, 
shear-shear correlations provide a powerful cross check 
on cosmological parameter determinations. It is well 
known that in principle shear-shear correlations are an 
extremely powerful probes of the dark energy parame- 
ters (e.g. |M|32). 

In Fig. |3 we show that with shear-shear correlations 
alone all three dark energy parameters w p i V ot, w a and 
Ode can be simultaneously measured. Note that we do 
not place external prior constraints on Ode unlike the 
convention commonly found in the literature but that 
the effect of a prior can be readily read off these figures 



For comparison with other probes, it is useful to note 
that the dark energy pivot point of the combined power 
spectra information is £ p ivot = 0.36 ~ ^de, the epoch 
of dark energy domination. Equation (|Allfl then maps 
the errors back into a more conventional description 
such as w (a n — 1). For example if we add in pro- 
jected CMB constraints for the Planck satellite [U, the 
errors improve to o-(SIde) = 0.004, a(w p i vot ) = 0.01 
[<j(w ) = 0.04], cr(w a ) = 0.08 as shown in Fig. H Note 
that the combined pivot point shifts to higher redshift. 
Here we have employed the sensitivities in [33| amount- 
ing to e.g. cr(lnD*) = 0.002 where is the angular 
diameter distance to last scattering; note that current 
constraints are at the 0.04 level. 

Finally, with all three power spectra one can probe the 
evolution of the halo parameters and hence aspects of the 
evolution of the underlying galaxy population. We show 
in Fig. [SI the resulting errors on the 5 halo parameters 
as a function of redshift employing all 2-point informa- 
tion out to lp — 3000. The errors for different redshifts 
are nearly independent whereas the errors between the 
5 halo parameters at a given redshift are strongly cor- 
related. The correlation indicates degeneracies between 
the parameters as can also be seen in Fig.^J Thus there 
are combinations of halo parameters that are better de- 
termined than implied by Fig. [SJ These combinations 
control the shape and amplitude of the spectra. For ex- 
ample, the linear combination that controls the bias and 
correlation coefficient across the 0.3 < z < 0.4 bin and at 
scales central to the constraint, e.g. k = 0.3h Mpc -1 are 
separately constrained at the level of cr(ln6) = 0.01 and 
a(\nR) = 0.004 (fiducial values: b = 2.2 and R = 1.2). 
Constraints on the correlation coefficient mainly reflect 
its insensitivity to halo parameters on large scales; they 
weaken in the deeply non-linear regime. Constraints on 
the bias as a function of scale remain at the percent level 
from the linear to well into the non-linear regime. 
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FIG. 5: Halo parameter errors as a function of redshift for the 
Nl = 10 lens galaxy bins and Ns = 5 source galaxy bins with 
all 2-point information to lp — 3000 and the efficiency ratio 
on all scales. Parameters at different redshifts are largely 
uncorrelated whereas those at the same redshift are highly 
correlated. Linear combinations that control the shape and 
amplitude of the power spectra are better constrained (see 
text). 



IV. DISCUSSION 

The joint analysis of galaxy clustering and lensing data 
from next-generation surveys offers unique opportuni- 
ties to simultaneously determine the evolution of clus- 
tering in the dark matter and galaxies. These surveys 
are expected to provide multi-color catalogs of galaxies 
with well characterized photometric redshifts well beyond 
z = 1. We have shown here that with redshift informa- 
tion, even at the 2-point or power spectrum level, there 
is enough information in the galaxy-shear correlation, es- 
pecially when combined with galaxy-galaxy correlations, 
to jointly solve for a model with ~ 50 parameters to 
describe the galaxy evolution and 3 dark energy param- 
eters. We have conservatively allowed the galaxy param- 
eters to vary independently as a function of lens redshift 
and marginalized the associated parameters when quot- 
ing dark energy constraints. 

The dark energy parameter determinations are statis- 
tically competitive with the shear-shear correlations and 
should be more robust to systematic errors in the shear 
determinations. They furthermore can provide a better 
redshift localization of dark energy effects given the broad 
lensing kernel of the shear. 

These determinations are assisted by two relatively ro- 
bust features in the halo model: on large scales, the com- 
bination of galaxy-mass correlations and galaxy-galaxy 
correlations may be used to determine the underlying 
mass-mass correlations in a manner that is only weakly 
sensitive to halo model assumptions; on small scales 
the ratio of galaxy-shear correlations at different source 
redshifts yields distance ratio information that is only 



weakly dependent on the evolution of the galaxy-mass 
correlation. 

Equally important, by combining galaxy-galaxy and 
galaxy-shear clustering one can determine whether the 
halo model employed here suffices as a description of 
the relationship between galaxies and dark matter ha- 
los. Consistency of the halo determinations can also be 
checked by selecting lens galaxies with multiple cuts on 
luminosity or rarity. Consistency of the dark energy de- 
terminations can be checked against the shear-shear cor- 
relations which depend only on the mass spectrum. A 
further extension to include contributions from the many 
galaxy-shear bispectra would also improve parameter ac- 
curacies and provide cross-checks [34|, |35j . 

Halo model parameter determinations will be valuable 
in testing models of galaxy formation (e.g. |3o. 1371 I3q1. 
Galaxy parameters measured from surveys as a function 
of galaxy type, luminosity and redshift can be compared 
with simulations and semi-analytic model predictions of 
the same parameters. Although these studies may show 
that the our halo model requires further modification and 
extension, we believe that the prospects are bright for a 
joint solution of galaxy and dark matter clustering. 
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APPENDIX A: COSMOLOGICAL POWER 
SPECTRUM 

The linear mass power spectrum P(k) which under- 
lies all observable power spectra depends only on well- 
motivated cosmological parameters in the context of the 
successful ACDM cosmology. In this section, we describe 
its parameterization in terms of the initial conditions HA II 
and evolution HA 21 

1. Initial Conditions 

We begin by assuming that massive neutrinos make 
a negligible contribution to the matter density. The 
shape of the mass power spectrum is then specified by 
the baryon density flbh 2 , the dark matter density fl m h 2 , 
and the spectrum of initial curvature fluctuations £ 

A ? = < A1 > 

where fco = 0.05 Mpc -1 is the normalization scale. We 
take fiducial values for these parameters that are con- 
sistent with the WMAP determinations: VL b h 2 = 0.024, 
fl m h 2 = 0.14, n = 1 and S c = 5.07 x 10~ 5 52]. The 
current uncertainties in these parameters are at the 10% 



10 



level or better. Note that 5q is related to the WMAP 
normalization parameter by 



4\2 



A = (1.84(5 C x 10 4 ) 



(A2) 



and current and future uncertainties in this parameter 
are expected to be dominated by uncertainties in the 
Thomson optical depth to reionization r, i.e. 



5 C m 5.07e- (ai7 - T) x 10~ 5 . 



(A3) 



Thus the power spectrum as a function of k in Mpc -1 
(not h Mpc -1 ) in the matter dominated regime can be 
considered as largely known. 



2. Evolution 

The shape of this initial power spectrum does not 
change during dark energy domination on scales below 
the sound horizon of the dark energy. The amplitude of 
the linear power spectrum depends on the initial normal- 
ization 6$ and the "growth function" , here the decay rate 
of potential fluctuations G 



P{k,z) 



1 G(z) 



1 + z G 



P(k,0). 



(A4) 



where Go = G(z = 0) and we assume that all relevant 
scales are sufficiently below the maximal sound horizon 
of the dark energy. P{k, 0) can be evaluated from any 
one of a number of Einstein-Boltzmann codes. 

The normalization of the linear power spectrum today 
is conventionally given at a scale of r = 8/i _1 Mpc 

d 3 k 



{2 ^P(k,0)W 2 (kr) 



5.59 x 10~ 5 V 0.024 



-1/3 



0.14 



x (3.123/i) 



(n-l)/2 



0.72 



0.693 



Gp 

0J6 ' 



(A5) 



where W a [x) — 3a; -3 (sin x — xcosx) is the Fourier 
transform of a top hat window. The approximation in 
Eqn. l)A5j) is valid to the 1% level for individual variations 
of the parameters in the regime 0.019 < f2&/i 2 < 0.03, 
0.11 < Q m h 2 < 0.18, 0.7 < n < 1.3, 0.5 < h < 1 which 
more than span the current observational errors. Note 
that because the normalization is given in hr x Mpc, there 
is a strong scaling with the Hubble constant. This scal- 
ing actually assists dark energy determinations since in a 
flat universe IIde = 1 — and precise measurements of 
fi m /i 2 makes h depend on Qde only. A measurement of 
the Hubble constant is a measurement of the dark energy 
density. Likewise a measurement of o$ is a measurement 
of a specific combination of dark energy parameters. 

We will hereafter limit ourselves to flat universes but 
correspondingly neglect the dark energy information from 



the CMB unless otherwise specified. The angular diam- 
eter distance to recombination D* has been measured 
to aQnD*) = 0.04 and this constraint will continue 
to improve with better measurements of Q m h from the 
peaks as (see e.g. 0) 



cr(lnD*) w -a{\nn m h 2 ) . 



(A6) 



The rationale behind dropping this constraint is that this 
measurement will be used in conjunction with galaxy and 
lensing constraints to eliminate any small curvature con- 
tribution that might exist. In a flat universe, closely 
follows Go in its dark energy dependence and so may be 
used as a powerful consistency test for the absence of 
spatial curvature. We discuss this issue further in TOI CI 
The growth function G then depends only on the dark 
energy density Qde(cl) = 8itGpde/3H 2 and equation of 
state w(a) = pde/pde through the equation (e.g. |33l |') 



d 2 G 
din a 2 



dG 
din a 



+ -[1 - w{a)}n BE (a)G = . 



(A7) 



where the initial conditions are G = 1 and dG/d\na = 
at an epoch well before dark energy domination z 3> zde, 



1 + zbe = a 



DE 



DE 



1 - fi 



DE 



-l/3ioo 



(A8) 



where Qde = ^de(z = 0) and wq = w(z = 0). 

Given a constant equation of state, Go follows the ap- 
proximate form 



Go 
F(x) 



0.76 



- - 1< 



0.236 



v 0.27 
(1 + 0.498a; 



F[n 4 ^(l + wbe) 



(A9) 



where the approximation holds to ~ 1% for separate vari- 
ations of —1 < wde < —1/2 and 0.12 < f2 m < 0.5. Here 
and throughout we denote a dark energy parameteriza- 
tion for which w is constant with u>de = w. Note the 
fairly strong scaling of Go with wde around the fiducial 
model AGo/Go « -0.33Aw D e. 

A dynamical form of dark energy is unlikely to possess 
a strictly constant equation of state. Since it only has 
observable consequences at z < zde, it is convenient to 
describe the function with its first order Taylor expan- 
sion. We therefore choose an equation of state parame- 
terization 



w(a) = w n + (a n - a)w a , 



(A10) 
this general- 



with the expansion around some epoch a v 
izes the form employed in |4lj | where a n = 1. 

We instead choose a n = a p i vo t such that the errors 
in Wpivot and w a for a given observable are uncorrelated 
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•"-"pivot 



0. The pivot point a p iv t can be derived from 



10 



a general representation via the transformation 



OC0 



dpu 



(All) 



For the transformation to the pivot representation [4: 



dw„ 



1, 



dw a 



^pivot ) 



(A12) 



from which it follows that the errors decorrelate for a 
shift in the pivot of 



C-u. 



*>n ^pivot 



(A13) 



Moreover the resulting errors on w p i vo t are then equal 
to those on wde, a constant w. The pivot redshift is 
therefore also the redshift at which w is best constrained. 

The drawback to choosing the pivot redshift is that it is 
specific to the observable, survey and dark energy model. 
As the redshift where the dark energy evolution is best 
constrained, the pivot point for growth measurements 
tends to coincide roughly with z n ss 0.4 or zde in the 
fiducial A model. In fact, further taking u>de — > w(sde) 
in Eqn. (|A9|) yields an approximation to the growth func- 
tion that is good to several percent across a wide range 
of w a . Choosing this standard normalization epoch then 
provides the benefits of the pivot point without the draw- 
backs. Although we will employ the pivot redshift for the 
galaxy-lensing study below, it is sufficiently close to zde 
that it may be interpreted in this manner. Note that the 
pivot point for distance-based dark energy measurements 
can be at even higher redshift zde ^ z < 1 so that it be- 
comes even more important to choose a z„ ^ for the 
characterization of constraints. 

In summary, our linear matter power spectrum is spec- 
ified by 7 parameters: 4 that are already well constrained 
by the CMB £ f (= 5.07 x lO" 5 ), n{= 1), n b h 2 (= 0.024), 
fl m h 2 (— 0.14), and 3 dark energy parameters fi DE (= 



0.73), 



J pivot ^ 



T) (or u>q), w a (= 0) which galaxy- 



lensing power spectra can help constrain. Parameter val- 
ues of the fiducial cosmology are given in parentheses. 
Given a linear power spectrum and cosmology, cosmo- 
logical simulations can accurately predict the fully non- 
linear mass power spectrum and hence the shear-shear 2- 
point correlations. On the other hand, the power spectra 
involving galaxies will require some semi-analytic mod- 
eling for the foreseeable future. We now turn to a halo 
model for the parameterization of the underlying rela- 
tionship between the observable power spectra and the 
linear theory predictions. 



APPENDIX B: HALO MODEL 
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FIG. 6: Host halo mass function dn^/dlnM^ and halo occupa- 
tion distribution function N g (Mf l ). Three halo parameters control 
the shape of the distribution: A B the satellite-host normalization 
or crossing point, m s the satellite slope (which pivots N a about 
the crossing point) , and o\ n m the scatter in the mass observable 
relation. Galaxies selected by a flux limit are matched in num- 
ber density by adjusting the threshold mass iWthi here illustrated 
for M t h = 10 13,5 /i — 1 Mq. Here and throughout the fiducial model 
is a flat ACDM model with cosmological parameters f2 m = 0.27, 



h = 0.72, n s = 1, St = 5.07 x 10" 



and halo parameters A B 



- 30, 



(cr 8 = 0.91) and Q b h 2 = 0.024 
= 1 and (T\ n a/ = 0.1. 



underlying three dimensional galaxy and mass density 
power spectra. Recent work on comparing simulations to 
galaxy clustering data have shown that to first approx- 
imation, galaxies selected on luminosity are assigned to 
a mass-based selection of dark matter (sub)halos of the 
same spatial number density [l2T |. This mapping avoids 
the traditional problem of defining explicit halo mass- 
luminosity relationships. In ^F3 II we build this underly- 
ing ansatz into a halo description of the galaxies and mass 
(see |43| for a review) and obtain a parameterized model 
for their joint power spectra in %B 21 We then describe 
the phenomenology of the fiducial halo model in 3l and 
study the sensitivity of the power spectra to variations 
in the halo model and cosmological parameters in %B 41 



1. Host and Satellites 

Under the halo model, the power spectra are described 
by the abundance, clustering, density profile, substruc- 
ture and galaxy occupation of dark matter halos. For 
the comoving abundance of host halos, we take the mass 
function E3 



To study the information contained in the lensing and 
galaxy two point observables, we require a model for the 



dn h _ po dv 
dlnM ~ M J[V ' d\nM 



(Bl) 
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where po is the matter density today po = p m (z = 0), 
v = 5J<7 



f(v) = A^av 2 [l + [av 2 )- p ] 



cxp[-at/ 2 /2] ■ (B2) 



Here cr(M; z) is the rms of the linear density field 
smoothed with a top hat of a radius that encloses the 
mass M . We choose S c — 1.69, a — 0.75, p — 0.3, and A 
such that J dvf(v) = 1. The mass function of the fidu- 
cial model is shown in Fig. bracketing the redshifts of 
interest for lensing (0 < z < 1). 

The halo cluster ing is then given by the peak- 
background split as [13 



b(M) = 1 



av* - 1 



2p 



(B3) 



S c 6 c [l + (av 2 )P}' 
For the halo density profile, we take the NFW form ^1 

( M \ _ Avp0 - 1 

P[r ' ,C) 3 ln(l + c) - c/(l + c) ifc(l + i?c) 2 ' 

(B4) 

where R — r/r v with the virial radius 



= / 3M„ 

\47I-jOo 

We take the concentration |H| 



c(M v ) = -L- ( ^ 



1/3 



-0.13 



(B5) 



(B6) 



where M* is defined by er(M»; z = 0) = <5 C . We convert 
between the halo mass M, assumed to be defined at an 
overdensity of 180 times the mean density, and the virial 
mass M v defined at an overdensity |49j 

18^ 2 + 82» DE (z)-39[» DE (z)] 2 , R7 ^ 
AZ) ~ 1 + SlMz) ' [ ' 

using the NFW profile (see e.g. [HoT]). 

The profile enters into the power spectra through its 
normalized Fourier transform 



y{k,M, c) 



1 sm ikft 

— dr47rr 2 p(r 7 M 7 c)^-^. (B8) 



M 



kr 



Although there are current uncertainties in these descrip- 
tions of the halo mass function, bias and concentration, 
we choose not to associate free parameters with these 
dark matter clustering properties. The characterization 
of these properties and more importantly the mass power 
spectrum itself can in principle be fixed by better simu- 
lations. Similarly for the association with galaxies, since 
we will be matching number densities of objects, the halo 
mass definition employed here may be replaced with any 
variable that defines the selection of objects in the simu- 
lations and is a monotonic function of galaxy luminosity 
on average. 
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FIG. 7: Satellite fraction dris/dfoh + n s ). The satellite fraction 
at a given mass for a given halo occupation distribution N(Mfo) is 
largely a function of the peak height threshold v in Eqn. IB1I as 
rare objects are unlikely to be satellites. The conversion to mass 
at z = is given in the upper axis. The range of fractions for 
M = 10 13 - 5 h- 1 M e and < z < 1 is shown shaded. 



The halo model predicts the galaxy and galaxy-mass 
power spectra under an assumption for the statistics of 
the occupation of the host halo by galaxies [H ■ Since 
each galaxy also carries its own dark matter halo we will 
hereafter distinguish between the host halo of mass Mh 
and its satellite halos of mass M s . For simplicity we 
take the satellites to also have NFW profiles but with 
an adjustable concentration; a more sophisticated model 
would account for the change in the functional form due 
to truncation from tidal stripping and trends in mass. 

We begin by taking a simulation-motivated universal 
form for the number of dark matter satellites in the host 
halo as a function of their mass ratio |53l . IH^j 



dN a 
din M s 



Mh 
AM, 



(B9) 



for Mh/ (100 A s ) < M s < Mh where simulations suggest 
that A s w 30 and m s « 1 [l3] • The cutoff to low masses 
prevents a mild logarithmic divergence in the total mass; 
this arbitrary cutoff only affects the mass power spectrum 
at unobservably small scales. 

We next associate galaxies with these satellites. The 
total number of satellite galaxies above a given threshold 
mass M t h then becomes 



dM s dN s 
Mt h M s d\nM s 
M h 



A S M, 



ih 



(B10) 



To this population of satellite galaxies we add the central 
galaxy associated with the host halo itself to obtain the 
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FIG. 8: Relative galaxy-galaxy and galaxy-mass power spec- 
tra. In the linear regime, k < O.lh Mpc -1 , both b = 
(pS 9 s g / p S m 6 m y/2 and b / R = p s g s g i p 5 g s m return the lin _ 

ear bias which increases as the objects become rarer. In the 
non-linear regime, P s s s a anc j p s g s m con tinue to track each 
other and exceed the mass power spectrum. 



total number of galaxies above the threshold (see Fig. El 



N g (M h ;M th )=N h +N a 



(Bll) 



These two pieces have different statistical properties. 
The central galaxy may be either above or below thresh- 
old and hence either occupied or unoccupied leading to 
(N%) = (Nh) — Nh- We model the average value as a 
step function at some limiting mass M tn smoothed by a 
Gaussian in lnM^ to reflect scatter in the conversion of 
a magnitude limit to a mass limit 



N h = ^Erfc 



lnM th /M h 
V2< 



O In M 



(B12) 



The dark matter satellites follow a Poisson distribution 
with (Nf, = N 2 a + N B 02. 

Given a fixed shape for N g , the threshold mass M t h is 
not a free parameter. Rather it is fixed by matching the 
space density of galaxies to the space density inferred by 
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FIG. 9: Correspondence of spatial wavenumber k (h —1 Mpc) to 
angular wavenumber I as a function of redshift z given by the Lim- 
ber approximation in the fiducial cosmology. Also shown is the 



projected non-linear scale k n i where k^P 



1. Note that 



I ~ 10 corresponds to the mildly non-linear regime for the red- 
shifts in question. 



the observed number counts 



n v (z;M t h) = 



dMh , r , , , dn h 



Jo Mh 

n A W g (z) 



■N(M h ;M th ] 



dlnM h 



(B13) 



Note that by choosing rare objects, the implied M tn > 
M* and so the population will be dominated by host 
galaxies rather than satellite galaxies (see Fig. EJ). This 
fact will be useful in minimizing the uncertainties and 
systematic errors associated with the model for N s and 
the satellite profiles, e.g. the mass scaling of the con- 
centration parameters. The satellite contribution can be 
quantified by considering the satellite mass function 



dn s 



din M s 



dM h dN» 



dnh 



M, 



Mh din M a dlnMh 



dM h 
M h 



M, 



h 



A S M S 



dnh 



din Mi 
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and comparing it to the host halo mass function itself 



dn x 



dn„ 



d{rth + n s ) din M 



d{n s + n h ) 
dlnM 



(B15) 



This ratio depends mainly on the peak height threshold 
v for a given form for the shape of N g (Mh) (see Fig.|7J. 
With our fiducial choice, it saturates at the low mass end 
at ~ 30%. With a selection of rare objects, e.g. corre- 
sponding to Mth ~ 1O 13 5 /i _1 M0 in the fiducial model, 
the satellite fraction can be less than 10% at all redshifts. 
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2. Power Spectra 

The power spectra are defined by the host and satellite 
distributions as 



jS„S 



(k) = I 2 lm (k)p(k) + I 2m (k) , 

r d °(k) = i 2 lg (k)p(k) + i 2g {k) , 

(k) = h g {k)h m {k)P{k) + h c {k) . (B16) 



The indices 1 and 2 refer to the number of points in a 
single halo. The first piece then represents two points in 
separate host halos correlated by the linear power spec- 
trum P{k) 



dM h ( M h 
M h \ po 
J_ f dM h 

fly 



dni-, 



dhxM, 



[N h + N s y g ] 



-b(M h )y h , (B17) 



dn h 
d\nM h 



b(M h ) , 



and the second piece is the contribution from two points 
within a parent halo including its satellite contributions 



h g 

he 



dM /MY 
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dn 
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2N h N s y g ] 
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d\nM 

dn h 



d\nM 



dm* 



d\nM 



H(M - M th )y s 



dlnM ' 
{NhVh + N s y h y g ) 
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Here we have employed the shorthand convention j/^ jS = 
y(k,Ch, s , M h , s ), V g = y(k,c g ,M h ) and the step function 
H(x) = 1 for x > and H(x) = for x < 0. The effect of 
satellite galaxies occupying satellite halos takes the same 
form as the central galaxy occupying the host halo 01 
except that we have neglected the small effect of smooth- 
ing around the threshold mass. In this description, the 
satellite-satellite mass correlation and satellite-halo mass 
correlation are implicitly included as the replacement of 
mass lost to subhalos in the term involving the parent 
profiles yh in Eqn. (|B18|) . 




FIG. 10: Angular power spectra in the fiducial halo and 
cosmological model for two galaxy (gi : < z < 0.5; 
g 2 : 0.5 < z < 1; M t h = 1O 13 ' 5 /i _1 M ) and source bins 
(ei : z < 1.5; €2 : z > 1.5; z mc d ; s = 1-5). Note the inflec- 
tion at I ~ few xlO 2 and the large number of cross spectra 
(N L N S = 4). 



0.01 
0.1 



R (h l Mpc) 

0.1 1 



A 
& 



0.01 



0.001 







\ ^^0 14 /r'M o ; 




^\10 13 \. 




\io 12 \ \ _ 


■z L =0.5;z s =1.5 

1 . . , . 


1111 ■ 



100 



CI 

o 
Q. 



10 t © 



< 



0.1 1 
9 (arcmin) 



10 



Fiducial Model 



FIG. 11: Tangential shear and surface mass density differ- 
ence in the fiducial halo and cosmological model for several 
choices of the mass threshold. Here we take a delta function 
distribution of lenses (zl = 0.5) and sources (zs = 1.5). Note 
the turnover on small scales (0.5 — 1') and inflection near 10' 
causing excess signal at large angles. 



Five functions of redshift specify our halo model, the 
satellite- host normalization A s , the slope of the satel- 
lite mass function and occupation distribution m s , the 
scatter in the mass observable relation a\ n Mi the con- 
centration of the satellites in the host halo c g and the 
concentration of the mass profiles of the satellites c s . We 
assume that the latter two follow the mass scaling of iso- 
lated halos in Eqn. 1|B6(1 but have an arbitrary normal- 
ization. We further represent the functions by a set of 
5Nl parameters that specify their values at the redshifts 
of the Nl lens redshift bins. 



The values of these parameters will in the future be de- 
termined by fitting to the joint power spectra. However, 
to study the potential of these future data sets we must 
specify their values in the fiducial model. The sensitivity 
of observables to variations in their values around this 
model is then quantified through the Fisher matrices of 

arm 

We choose fiducial values that are roughly consistent 
with low redshift measurements from the SDSS lensing 
data and TV-body simulations: A s = 30, m s = 1, 
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cinAf = 0.1 for all redshift bins. These three parameters 
define the shape of N g shown in Fig. [S] In our fiducial 
model the concentration parameters c s — c g = c^. 

In Fig. |S] we show the galaxy-galaxy and galaxy-mass 
power spectra with the fiducial halo and cosmological 
parameters. They are shown relative to the mass-mass 
spectrum at z — for halo abundances corresponding to 
several different mass thresholds in the fiducial ACDM 
cosmology (see flXjl . In the linear regime k < O.lh Mpc" 1 
the ratios 



and 



b(k) = 



b(k) 



pS g 5g 

pS m 5 m 



pSgSg 



1/2 



R(k) P s 3 



(B19) 



(B20) 



both return the constant linear bias of the objects. In 
other words, in the linear regime the correlation coeffi- 
cient 



R(k) 



8 a 8„ 



rps g s B ps m s m y/2 



(B21) 



independent of the population or halo model parameters. 
In this regime, a measurement of P s n Sm and P s s s s is a 
measurement of the mass power spectrum P and hence 
in combination can be used to study the dark-energy de- 
pendent growth rate. The linear bias increases with the 
mass threshold M t h or more generally with the rarity of 
the objects. 

In the non-linear regime, the bias becomes strongly 
scale dependent reflecting the correlations within a par- 
ent halo. On the other hand, the ratio b/R remains 
remarkably constant Formally the combination in- 
dicates that the galaxy-mass correlation R > 1. Un- 
der the halo model, every galaxy has a dark matter halo 
around it and so the relevant power spectrum for comput- 
ing the correlation coefficient is corrected for the excess 
shot noise power p s a s a + fiy 1 . The scale dependence of 
b and b/R can be used to pin down the halo parameters. 
For fixed halo parameters, it marks the non-linear scale 
which also depends on the dark energy through the linear 
growth rate. 

These properties remain qualitatively true for angular 
power spectra with the caveat that projection effects can 
broaden the non-linear transition regime for wide red- 
shift shells for the lens galaxies and the broad efficiency 
for lensing. In Fig. [5] we show the correspondence of k 
and I in the fiducial cosmology and under the Limber 
approximation. In Fig. I1UI we show an example of the 
Limber projection to the angular power spectra. Note 
that the inflection caused by the transition between the 
one host halo and (linear) two host halo regimes occurs 
at wavenumbers of I ~ few x 10 2 . Note that there are 
Nl = 2 galaxy spectra, N^Ns = 4 cross spectra, and 
Ns(N s + l)/2 = 3 shear spectra. 
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FIG. 12: Halo parameter sensitivity compared with dark energy 
sensitivity of galaxy-galaxy and galaxy-mass power spectra. In 
the non-linear regime, galaxy-galaxy power spectra become highly 
sensitive to most halo parameters. Galaxy-galaxy power spectra 
are insensitive to dark energy effects on the growth rate since a 
lowering of the amplitude is compensated by an increase in bias 
due to the increased rarity of the objects. Conversely the galaxy- 
mass power spectra are less sensitive to halo parameters and more 
sensitive to the dark energy. 



In Fig.EI we highlight the one halo regime by showing 
the predictions for the tangential shear from Eqn. 115|) 
for a given source and lens configuration. Note that the 
shear turns over as the angular radius resolves the scale 
radius of the halo and shows another break on large scales 
marking the beginning of the two halo regime. 



4. Parameter Sensitivity 

The sensitivity of the power spectra to halo and cosmo- 
logical parameters is quantified by their derivatives with 
respect to the parameters in Eqn. (|25() . In Fig. ^] we 
compare the halo parameter sensitivity of galaxy-galaxy 
and galaxy-mass spectra to the equation of state param- 
eter wde at z = for a population with an abundance in 
the fiducial model corresponding to M t h = 10 5 h~ 1 M Q . 
Since the initial normalization of the power spectrum is 
here fixed, sensitivity to iude is equivalent to sensitivity 
to the present-day normalization as (see Eqn. (|A9() 1 . 



16 



Variable Definition 


Eq. 


N L 


Lens galaxy bins 




N s 


Source galaxy bins 




N h 


Host halo occupation 


1BH1 


N s 


Satellite halo occupation 


ijmoji 


N a 


Total halo occupation 


1BH1 


M 


Halo mass 


m 


M h 


Host halo mass 




M s 


Satellite halo mass 


(lB9l 


M v 


Virial mass 




ri 


Angular position 


© 


n 


Initial tilt 


ED 


riA 


Galaxy angular density 


© 


riv 


Galaxy space density 


GJ 


rih 


Host halo space density 


(TeTTI 


n s 


Satellite halo space density 


ilBTH 


w 


Dark energy equation of state 




WDE 


Constant w(a) 




Wn 


Specific epoch w(a n ) 


JXToTi 


™pivot 


Best constrained w(a p i vo t) 
Present day w(a = 1) 




Wo 


lEl 


W a 


Evolution —dw/da 





TABLE I: Easily confused variables. 



The galaxy-galaxy spectrum is strikingly insensitive to 
wde compared with the halo parameters that control 
Ng(Mh): A Sl m s and ai n A- Changes in the parame- 
ters that control its shape change the power spectrum 
on all scales since they enter into the calculation of the 
mass threshold Af t h and hence the bias of the host ha- 
los. A change in the cosmology that for example low- 
ers the present amplitude of the matter power spectrum 
makes galaxies of a given number density rarer. Since 
rarer objects are more highly biased tracers of the mat- 
ter, the galaxy-galaxy power spectrum remains largely 
unchanged. The same is not true for the galaxy-mass 



power spectra. Indeed the galaxy-mass power spectra are 
substantially more sensitive to cosmological parameters 
than halo parameters. The combination of galaxy-galaxy 
and galaxy-mass power spectra is then particularly pow- 
erful for simultaneously determining the halo and cosmo- 
logical parameters. 

Neither set of spectra are very sensitive to the concen- 
tration parameters c s and c g except at very k values that 
correspond to I » 10 3 for the redshifts in question (see 
Fig. 0) ■ By definition, the concentration parameters only 
affect the power spectrum on small scales. Furthermore 
since we have chosen a high mass threshold, the fraction 
of objects that are satellites and hence affected by these 
parameters is low. This insensitivity helps justify our 
crude treatment of the profiles. 

As for the mass-mass power spectrum, and hence the 
shear-shear power spectrum, it is of course the most sen- 
sitive to the dark energy parameters being directly re- 
lated to the linear growth function of Eqn. i|A7() . It for- 
mally also depends on the satellite mass function and 
hence the halo model parameters that control it. How- 
ever we here take the perspective that the form of the 
mass-mass power spectrum as a function of cosmology 
will in the future be fixed directly by simulations re- 
placing the halo model description. Hence we take the 
uncertainties in the halo parameters to only affect the 
galaxy-galaxy and galaxy-mass power spectra. Opera- 
tionally we assume that variations in P" 5 ™" 5 ™ due to the 
variations in the satellite mass function are localized to 
the galaxy mass regime and compensated by variations 
at other masses to leave the spectrum invariant. This 
corresponds to dropping the parameter derivative terms 
in the Fisher matrix H24|) . 

We provide a guide to easily confused variables in 
Tab.Hl 
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